The interaction of light with an effective moving dispersive medium. 
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Intense laser pulses excite a nonlinear polarisation response that may create an effective flowing 
medium and, under appropriate conditions, a blocking horizon for light. Here we analyse in detail 
the properties and interaction of light with such laser-induced flowing media, fully accounting for 
the medium dispersion properties. A first Born approximation scattering model is found to be in 
excellent agreement with numerical simulations based on Maxwell's equations. Diamond is proposed 
as a promising candidate medium for future studies of Hawking emission from artificial, dispersive 
horizons. 



Recent developments in the understanding of light 
propagation have led to evidence that by using intense 
laser pulses propagating in a nonlinear medium, it is pos- 
sible to create an effective medium that flows with the 
same speed as the laser pulse, i.e. at speeds close or even 
higher than the speed of light in vacuum [1-9]. Indeed, 
in a medium with a third order (also called "Kerr" ) non- 
linear polarisation response the refractive index of the 
medium is given by n = n + n 2 I(z — vt), where n is 
the background index, ni is the nonlinear Kerr index 
and I(z — vt) is the laser pulse intensity profile, travel- 
ling along the z direction with velocity v. In a typical 
condensed medium, e.g. glass, the maximum amplitude 
of the laser pulse induced refractive index perturbation 
is Sn 

max — ^ 

0.01 - 0.001. This is sufficient to 
scatter light, e.g. either light from the laser pulse it- 
self (self-scattering process) or from a second weak probe 
pulse (induced scattering process). As discussed in the 
following, there is also interest in the fact that under ap- 
propriate conditions one may potentially observe induced 
scattering (and excitation) of photons that originate from 
the vacuum state. 

One motivation for the interest in the interaction be- 
tween a moving refractive index variation, Sn, and the 
vacuum state fluctuations lies in the prediction that the 
Sn may be described in terms of a horizon that mimics 
the event horizon of gravitational black hole [1-7]. In- 
deed, in the absence of material dispersion, it is possible 
to write the effective spacetime metric for the moving Sn 
as ds 2 = c 2 dt 2 + (dr — Vdt) 2 . This is equivalent to the 
Painlcvc-Gullstrand metric showing that in the reference 
frame of the Sn the medium is effectively flowing with 
speed V, that is a Galilean velocity (—00 < V < +00) 
and is determined by n and v. A horizon is formed when 
V = c [5, 10]. Interaction with the vacuum state leads to 
the emission of photon pairs that are analogous to Hawk- 
ing emission from a black hole. Recent measurements 
claimed the observation of a spontaneous emission from 
a laser pulse induced moving Sn that appeared to have 
some of the features predicted for the analogue Hawking 
emission [11, 12]. However, these measurements have not 



been considered as conclusive and further evidence is re- 
quired [13]. 

In this Letter we investigate in detail the interaction be- 
tween a moving Sn and a probe light pulse. Our ana- 
lytical and numerical models fully account, without any 
approximations, for material dispersion, the full shape 
of the Sn and allow a comprehensive analysis of the in- 
teraction dynamics for varying Sn speed. These studies 
highlight new and unexpected interaction dynamics be- 
tween light and a moving or flowing medium: the input 
probe pulse is scattered into two output modes one with 
positive and the other with negative comoving frequen- 
cies [14]. In agreement with previous studies this leads to 
amplified emission [15] with the novelty that the noise of 
this amplifier has a constant temperature for all emitted 
modes (frequencies) thus resulting in a blackbody emis- 
sion when seeded from the vacuum state. The blackbody 
emission is predicted independently of the Sn speed al- 
though the total photon number increases dramatically 
close to the condition for which a complete blocking hori- 
zon is formed. It is in this regime that we propose a set- 
ting for future experiments in diamond that exhibits the 
required conditions to observe analogue Hawking emis- 
sion. 

Scattering from a moving Sn: A soliton propagating in 
a dispersive medium will shed light through a mecha- 
nism known as resonant (or dispersive wave) radiation 
(RR) [17-20]. This emission may be described as a self- 
induced scattering process whereby light from the soliton 
is scattered into a frequency-shifted mode by the self- 
induced Kerr Sn. It was recently shown that a second 
scattered mode also exists: this mode has negative fre- 
quency in the reference frame comoving with the soliton 
and has been called "negative-frequency resonant radi- 
ation" (NRR) [14]. Both these modes are found within 
the framework of a model that neglects all nonlinear ef- 
fects and simply considers the soliton Sn and models how 
light is scattered within the first Born-approximation 
[15]. The model therefore generalises RR and NRR gen- 
eration beyond soliton physics to include also systems in 
that do not support solitons (e.g. a 3-dimensional pulse 
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with transverse spatial dynamics or an intense pulse in 
the normal group velocity dispersion regime). The scat- 
tered amplitudes of the two modes are given by [15] 

S M = S^V 1 "*^^ " ^)(1) 
2c 2 k z (wrr) 

S(-nrr) = - fnnno ^ ^ R {uJNRR + w IN )(2) 

where i? is the Sn Fourier transform and wrr^nrr.in are 
the RR, NRR and input probe pulse frequencies. A fur- 
ther important result is the generalisation of the photon 
conservation Manley-Rowe relations to the case of a mov- 
ing medium [15, 21, 22]: \RR\ 2 - \NRR\ 2 = 1, where 
| • | 2 indicates a photon number normalised to the input 
photon number. This generalised Manley-Rowe relation 
implies that the scattering process is actually a photon 
number amplifier, \RR\ 2 + \NRR\ 2 > 1. This has been 
confirmed numerically by monitoring the total photon 
number of a weak seed pulse that interacts with a soliton 
pulse and is subsequently scattered into RR and NRR 
modes [15]. 

Amplification and gain spectrum: In the following we 
present a series of numerical simulations and a compar- 
ison with predictions based on the Born approximation 
equations (1) and (2) with which we estimate the am- 
plitudes of the RR and NRR waves generated by an in- 
put probe pulse interacting with a moving Sn, for vary- 
ing frequency of the input pulse. The numerical sim- 
ulations were carried out with various codes i.e. using 
both FDTD [23] and the UPPE [15, 24] model all un- 
der the same input conditions and both delivered the 
same results. In Fig. (1) we show examples obtained with 
the FDTD code with Sn = <5n max exp [-[(z - vt) 2 jr 2 \ m \ 
where 5n max = 0.01, r and m are chosen so that the Sn 
front has a 7 fs rise time. The dispersion is chosen so as 
to resemble that of diamond and is shown in the figure 
in (u' ',u)) frequency coordinates, where primed quanti- 
ties indicate that they refer to the comoving reference 
frame. This is a convenient representation: momentum 
conservation that must be satisfied during the scattering 
process, translates in the comoving frame to a frequency, 
u/ conservation. All modes and relative laboratory fre- 
quencies, wrj^nrr may therefore be determined from the 
intersections of the dispersion curve with a horizontal 
line that passes through the input comoving frequencies, 
±<4n [15]- 

Figure 1(a) shows the near field, i.e. temporal profile of 
an input 5 /jm wavelength probe beam interacting with 
the Sn in the comoving frame for the generic case in which 
the Sn speed, v = 1.9 X 10 8 m/s is too slow to actually 
completely block the impinging light for <5n max = 0.01. 
We refer to this as the "non-blocking" case: the input 
light will be transmitted through the Sn and will only 
be partly scattered into the RR and NRR modes that 
are reflected backwards in the comoving frame (but are 
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FIG. 1: Example of a numerically simulated interaction of a 
probe pulse with a moving Sn. (a), (b) and (c) show the tem- 
poral profile evolution, spectral evolution and relevant disper- 
sion curves, respectively, for the case in which the 5n moves 
too slowly to form a blocking horizon, (d), (e) and (f) show 
similar figures for a faster Sn such that a blocking horizon is 
formed. The Sn rising front is show as a grey shaded area in 
(a) and (d). All intensity plots are shown over 4 decades in 
logarithmic scale. Stars in (c) and (e) indicate the position of 
the cut-off frequency, u; max . 



travelling forwards in the laboratory frame). This partial 
conversion is also clear in the frequency spectrum evolu- 
tion in Fig. 1(b) that shows generated peaks that are in 
excellent agreement with the frequencies predicted from 
the dispersion relation in Fig. 1(c). 

Figures 1(c), (d) and (e) show similar results but now 
v = 2.07 x 10 8 m/s has been increased such that the 
input light pulse is slowed down upon interacting with 
the Sn rising front to the point that its group velocity 
becomes equal to v, i.e. the Sn presents a blocking hori- 
zon for light. In this case the input mode is completely 
converted to the RR and NRR modes. 
We now examine the ratio r = | NRR\ 2 /\RR\ 2 for vary- 
ing input frequency. In Fig. 2(a) we show a typical ex- 
ample in logarithmic scale (red points are results from 
numerical simulations with a non-blocking v = 1.9 X 10 8 
m/s, the solid line is a linear fit). The surprising feature 
here is that r has a clear exponential dependence over 
more than 8 decades. If we then combine this expnential 
dependence r = exp(— aw') with the generalised Manley- 
Rowe relations we find that \NRR\ 2 = l/[exp(aw') — 1]. 




1 2 3 4 5 
to' (rad/ps) 





10°' 

r 




io; 2 r 

r 


+IH 


10 -4 r 


? 1 


io" 6 r 


NR 


r 

10" 8 r 




10" 10 f 



1.75 1.85 1.95 2.05 
v (10 8 m/s) 



1000 




0.1 0.2 0.3 0.4 



FIG. 2: (a) Numerical simulation for v = 1.9 X 10 8 m/s of 
r = \NRR\ 2 /\RR\ 2 for varying input frequency (dots) and 
best fit with exponential function (solid line). The Born 
approximation calculation is hows as a dashed line . (b) 
Temperature [derived from graphs as in (a)] for varying Sn 
speed - simulations (solid line) and Born approximation model 
(dashed line). The shaded area indicates speeds for which a 
blocking horizon is formed, (c) Total normalised photon count 
for varying speed, (d) Comparison between numerically esti- 
mated emission temperature (dots) and theoretical Hawking 
emission temperature estimated from the Sn gradient (solid 
line) for varying maximum index change, <5n max . 



This implies that the NRR mode emission follows a black- 
body law with temperature given by T = h/(k&a). We 
note that the backbody dependence is maintained also in 
the laboratory reference frame, albeit with a temperature 
?lab = T comoving /[7(l - vn/c)] where 7 = l/y/l - (v/c) 2 
[5]. We note that we can also evaluate r directly from 
the first Born approximation relations, Eqs. (l)-(2). The 
same exponential dependence for r is found (shown as a 
dashed line in Fig. 1(a)) and r — > 1 for win — > 0. 
In Fig. 2(b) we show the laboratory temperatures ob- 
tained from similar curves as shown in Fig. 2(a) for vary- 
ing Sn speeds: the solid line is the result from numerical 
simulations, the dashed line is the result from the Born 
approximation equations. The two curves are in over- 
all good qualitative agreement with even a quantitative 
agreement in terms of the actual temperature values and 
behaviour as the speed crosses over from the non-blocking 
to the blocking horizon case (shaded region in the figure). 
We note a dip in the temperature at the transition from 
non-blocking to blocking. This is due to the sharp cut-off 
of the negative mode scattering efficiency when the in- 
put frequencies intersect the dispersion curve close to its 
maximum value, w„ lax as predicted in Refs. [4, 25, 26]. In 
more detail, in the blocking case only the Sn — disper- 
sion curve [see e.g. Fig. 1(f)] has an w^ ax and suppression 
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FIG. 3: Numerical simulation of the nonlinear propagation 
of an 800 nm pump pulse and a co-propagating 15 THz probe 
pulse (shown in logarithmic scale over 4 decades). Temporal 
profile evolution of the 800 nm pump pulse (a), and of the THz 
pulse, (b). The arrow indicates the propagation distance at 
which a shock front forms on the pump pulse. "RR" indicates 
emission from the THz pulse scattered from the pump shock 
front, (c) shows the pump pulse spectral evolution. The THz 
pulse spectrum (6 decades in logarithmic scale) in (d) shows 
a clear peak centred around the u 1 — point, corresponding 
to a 7 rad/fs lab frame frequency. 



of the NRR mode is not observed as we are probing at 
relatively low frequencies. As v is decreased, the disper- 
sion curve relative to the internal region of the Sn raises 
up and a presents an u' max that lies in the same region 
as the probing frequencies. In agreement with previous 
studies, NRR generation is then strongly suppressed and 
the temperature decreases sharply [4, 25, 26]. If v is 
decreased even further the internal dispersion curve be- 
comes nearly identical to the external curve [similar to 
Fig. 1(c)] and the two-mode RR-NRR scattering process 
takes over again with a consequent return of high tem- 
perature blackbody emission. 

We now note that although the blackbody emission tem- 
perature appears to be high across all values of v, the 
emission efficiency is actually exponentially suppressed 
as v decreases. In Fig. 2(c) we show the total photon 
number \NRR\ 2 + \RR\ 2 evaluated from the numerical 
simulations. As can be seen, this number varies by sev- 
eral orders of magnitude with varying v and is dramati- 
cally enhanced in the presence of a blocking horizon. 
In the presence of a blocking horizon we may also evalu- 
ate the Hawking temperature (at the horizon) according 
to T = (l/2ir)"/ 2 v\dn/dt\ [5]. We find that the temper- 
ature evaluated both from the numerics and the Born 
model are in very good agreement with the predicted 
Hawking temperature. This is shown in Fig. 2(d) where 
we compare the (comoving) Hawking temperature (solid 
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curve) and the numerically estimated temperatures (red 
dots) for Sn with a fixed v and varying <5rt max - As can 
be seen, the agreement is excellent over a wide range of 
values, including Sn maK that are experimentally accessi- 
ble through the nonlinear Kerr effect. 
It has already been pointed out that the moving Sn acts 
as an amplifier. From the results presented here it is clear 
that the amplification increases with the Sn speed and 
reaches maximum efficiency in the presence of a block- 
ing horizon. If the input probe pulse were to be reduced 
to the level of the quantum noise fluctuations, then the 
results shown above predict the spontaneous emission 
of a blackbody spectrum. In other words, the noise of 
the amplifier is characterised by a temperature that re- 
mains constant across all frequencies. This is therefore 
a very different kind of amplifier with respect to better 
known examples in optics, e.g. optical parametric am- 
plifiers that have a noise temperature that scales linearly 
with frequency, T = huj/kB [16]. Finally, we note that 
the blackbody dependence implies that the amplifier gain 
scales as 1/uj' for (J — > 0. In the comoving frame light 
is coupled between modes that have a constant or nearly 
constant lj' ~ 0. Looking at Fig. 1(f) we see that if the 
input mode has lab. frequency ui close to zero then the 
output mode will appear at high frequencies, typically 
in the UV region, in the laboratory frame. The ampli- 
fier therefore converts radiation from the low to the high 
(lab frame) frequency modes and the 1 /to' blackbody di- 
vergence will appear as a divergence (or peak) centred 
at the UV lab frame frequency. The high amplification 
gain when coupling between a/ ~ modes also indicates 
possible methods for efficiently observing these scatter- 
ing effects and possible applications. 
For example, in Fig. 3 we show numerical simulations of a 
weak 15 THz probe (20 urn wavelength) pulse interacting 
with a TkSapph laser pulse with 800 nm wavelength, 60 
fs duration and input intensity 80 TW/cm 2 , propagating 
in a 500 /im think diamond sample. These simulations 
were performed using two UPPE equations, one for the 
pump and one for the THz pulse that are coupled only 
through a nonlinear cross-phase modulation term in the 
THz pulse equation cx 2n 2 I (details of the code can be 
found e.g. in Ref. [15]). The self-phase modulation term 
cx n 2 I is included in the 800 nm pump equation but 
four wave mixing and third harmonic generation have 
been purposely neglected. There is no blocking horizon 
for the probe pulse but nevertheless the spectrum clearly 
shows relatively efficient scattering to the UV that occurs 
predominantly when a shock front, i.e. the steepest gra- 
dient forms on the pump pulse [indicated by an arrow if 
Fig. 3(a)] in agreement with Eqs. (1) and (2). Most inter- 
estingly, the output mode is peaked around the u/ ~ 0, 
corresponding to a lab frequency of 7 rad/fs (270 nm 
wavelength). In other words, laboratory reference frame 
frequencies in the THz or multi-THz region are already 
sufficiently low to excite the 1/u/ gain of the amplifier. 
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FIG. 4: Numerical simulation of a superluminal 5n in fused 
silica, (a) dispersion relations outside the 5n and inside at 
<5n max = 0.008. (b) r = \NRR\ 2 /\RR\ 2 for two different Sn 
rise times as indicated in the figure, r is peaked at the pre- 
dicted comoving zero frequency, indicated with a filled square 
in (a) . Dashed curve - result from Born approximation calcu- 
lation. 



This in turn provides indications of the spectral range 
of the noise fluctuations that may be spontaneously ex- 
cited by the Sn and lead to the spontaneous (i.e. seeded 
by thermal or vacuum noise) emission peak in the UV. 
At room temperature, the background blackbody pho- 
ton mode density at 15 THz is ~ 10 5 and thus signifi- 
cantly larger that the quantum vacuum noise mode den- 
sity equal to just 1/2 a photon per mode: spontaneous 
emission will be seeded by the thermal background. How- 
ever, it is sufficient to cool the diamond sample to ~ 30 
K to invert the situation such that thermal fluctuations 
arc dominated by more than 1 order of magnitude by 
quantum vacuum noise. 

As a final comment, we note that superluminal speeds 
may be obtained using Bessel pulses: in diamond if v > c 
the dispersion curves both inside and outside the Sn re- 
semble the curve for Sn max in Fig. 1(f) so there is only 
one intersection point a/ = cj[ n . Therefore only one mode 
(NRR) will be excited and there will be no amplification. 
However, in fused silica glass the dispersion relation ex- 
hibits a resonance in the infrared that allows coupling to 
two output RR and NRR modes or between two u>' = 
modes. The dispersion relations for fused silica relative 
to recent experiments using such Bessel-induced super- 
luminal Sn are shown in Fig. 4(a). If the probe input 
pulse starts from inside the Sn, then the pulse will start 
to lag behind, exit the Sn and, when traversing the Sn 
gradient it will excite RR and NRR output modes. If 
we repeat the procedure outlined above, i.e. we repeat 
the simulations for varying input wavelength and moni- 
tor r = \NRR\ 2 /\RR\ 2 , we obtain the results in Fig. 4(b) 
that shows curves for two different Sn rise times (0.6 fs 
and 2 fs). The dashed curve shows the result obtained 
for the 2 fs gradient pulse using the Born approximation 
and is in good agreement with the numerics. In all cases, 
we see that r is peaked at the ui' = frequency point 
on the Sn — dispersion curve [indicated with a solid 
square in (a)]. This indicates that spontaneous emis- 
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sion will have maximum gain at this "zero-frequency" 
mode, in agreement with measurements of spontaneous 
emission reported in Ref. [11]. Although a quantitative 
comparison with analogue Hawking emission is somewhat 
thwarted by the absence of a blocking horizon, which in 
turn poses a difficulty in making a direct comparison with 
a Hawking temperature [13], we see that the gain peak at 
uj' <~ in fused silica glass is essentially of the same na- 
ture of that observed in the blocking (and non-blocking) 
horizon case analysed in diamond. However, the latter 
material does seem more appropriate for future experi- 
ments aimed at observing and studying a faithful ana- 
logue of Hawking emission. 

Conclusions: a numerical and analytical evaluation of 
the interaction of a weak light pulse with a steep re- 
fractive perturbation leads to the emission of two output 
modes. This interaction may lead to an increase of the 
overall photon number, i.e. amplification of the input 
light pulse at the expense of the moving 5n. If the Sn 
is generated through the nonlinear Kerr effect by an in- 
tense pump pulse, then energy will be extracted from 
this pump pulse. The moving Sn may be tuned so as to 
form a blocking horizon for incoming light modes. In this 
regime, the process is greatly enhanced and in diamond- 
like media (i.e. media with no resonances at infrared 
wavelengths) the emission bears all the attributes of ana- 
logue Hawking radiation, namely a blackbody emission 
originating from the excitation of two coupled modes at 
a blocking, yet dispersive horizon. Optical horizons have 
been proposed for all-optical transitors [27] but other ap- 
plications of an effective moving medium to be considered 
in future work could be relatively efficient THz detection 
by frequency conversion in diamond (even without a hori- 
zon, as shown in Fig. 3) or generation of squeezed vacuum 
states in the UV region. 
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